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Abstract 

Abstract. We propose an approximate description of basic parameters (ra- 
dius, mass and oblateness) of general relativistic compact rotating objects in 
terms of the parameters of the static configuration and of the angular veloc- 
ity only. The representation in terms of static properties is derived using the 
condition of stationary equilibrium together with some phenomenological as- 
sumptions. The predicted radius and mass of rotating neutron star (described 
by some realistic equations of state) and strange star (described by the bag 
model equation of state) are compared with data obtained by numerical in- 
tegration of gravitational field equation. The obtained formulae also allow a 
simple derivation of the "empirical" equation relating the maximum rotation 
frequency ^l ma x of uniformly rotating star models to the mass and radius of 
the maximum allowable mass configuration of the non-rotating model. 
Pacs Numbers: 



*E-mail:hrspksc@hkucc.hku.hk 
^E-mail: tcharko@hkusua.hku.hk 



1 



I. INTRODUCTION 



Rotation is a basic physical property of astrophysical objects. Study of rotational prop- 
erties of stars can lead to some restrictions imposed on the nuclear equation of state of 
dense matter at densities larger than nuclear density. Oscillations of rapidly rotating stars 
can become unstable hence producing detectable gravitational wave emissions. 

In the past decades there were numerous attempts to construct analytic models for 
rotating perfect fluid bodies in general relativity. But obtaining an exact interior solution 
for a rotating body proved to be a formidable task.The first step in this direction was the 
fundamental result of Kerr [1] who obtained the solution for the vacuum domain, outside the 
rotating star. It took no less than three decades of investigations to obtain the first, highly 
idealized model of a general relativistic, thin rotating disk of dust, by Neugebauer and Meinel 
[2]. Various schemes have been developed for obtaining stationary and axisymmetric perfect 
fluid solutions of the gravitational equations like Petrov type D, local rotational symmetry, 
fluid kinematics, non-trivial Killing tensor, vanishing Simon tensor, electric magnetic Weyl 
curvature, lagrangian or static- stationary symmetry, geodesic eigenrays etc. (for a recent 
review of rotating perfect fluid models in general relativity see [3]). 

If in the case of rapidly rotating stellar configurations there are still many unsolved 
problems, a remarkable progress has been made in the study of slowly and rigidly rotating 
perfect fluid configurations. By casting the metric in the form 



ds ' = e » dt 2 _ e A rfr 2 _ r 2 



d9 2 + sin 2 9 (# - codtf] + O (fi 2 ) , (1) 



Hartle [4] obtained a formalism that proved to be very useful in many investigations of 
the rotational properties of stars. However, this model, which considers only first order 
corrections in Q, can not be used to compute models of rapidly rotating relativistic stars 
with sufficient accuracy. 

On the other hand there has been recently a considerable advance in the numerical 
understanding of rotating stars. Several high precision numerical codes are now avalaible 
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and it has been shown that they agree with each other to remarkable accuracy (see [5] for a 
review of recent developments in numerical study of rotatation, nonaxisymmetric oscillations 
and instabilities of general relativistic stars). 

The non-sphericity of rapidly rotating stationary stellar configurations and the compli- 
cated character of the interplay of the effects of rotation and of those of general relativity 
seem not to permit a simple universal description of rotating compact objects. However, 
Haensel and Zdunik [6] and Friedman, Ipser and Parker [7] have found a simple relation con- 
necting the maximum rotation frequency fl ma x with the maximum mass M^ x and radius 
Rmax °f the static configuration: 



with Cg a constant which does not depend on the equation of state of the dense matter. 
The value of the constant Cs has been obtained by fitting equation (|2|) with data obtained 
by numerically integrating the gravitational field equations. It is given by Cs = 7200s -1 [7], 
Cs = 7700s -1 [8] or by Cs = 7730s -1 [9]. The empirical relation (Q) has been checked for 
many realistic equations of state for neutron stars [8], [10]- [12] and for the case of strange 
stars where the empirical formula also holds with a very good precision, the relative devia- 
tions do not exceeding 2% [9]. Equation (0), obtained on the basis of analyzing numerical 
solutions of the gravitational field equations, provides an enormous simplification of the 
problem of dynamical effects of rotation because the solutions of the complicate general 
relativistic equations for a rotating star can be replaced with the solutions of the much sim- 
pler TOV equation. Some attempts to explain this empirical relation were not concluded 
with satisfactory result. Weber and Glendenning [13], [14] used numerical models of slowly 
rotating relativistic stars to show that the formula still hold, but with Cs = 7500s -1 . Also 
in the slow rotation limit Glendenning and Weber [15] derived a formula relating VL max to 
M max / Rfnax, in terms of the mass, equatorial radius, moment of inertia, angular momentum 
and quadrupole moment of the maximally rotating configuration only. But it is not clear 
how the formula @ follows from their results. Up to now a clear physical understanding of 
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this relation is still missing. 

On the other hand a universal relation between the maximum mass and radius of non- 
rotating neutron star configuration and the mass M™^ and radius R^lx °f the configuration 
rotating with VL max of the form 

jijrot /~i j, ir stat rjrot /~f jjstat /o\ 

lvl max ^M^maxi ^max ^ R^rnaxi V 3 ) 

has also been found [12]. In @ Cm and Cr are specific equation of state dependent con- 
stants, whose values have been calculated, for a broad set of realistic EOS, in [12]. Their 
mean values are Cm = 1-18 and Cr = 1.34 [12]. The empirical constant Cs can be obtained, 
within an approximation better than 5% , from the formula [12] 

The general validity of (0) and (|3]) suggests the possibility that all basic physical param- 
eters of general relativistic rotating stellar objects (like mass and radius) can be somehow 
related to the similar parameters of the static configuration. It is the purpose of the present 
paper to propose a general description of basic physical parameters (mass, radius and oblate- 
ness) of compact general relativistic objects in terms of the physical properties of the static 
configuration and of the angular velocity only. To obtain the representation in terms of 
static physical parameters we use only the general relativistic conditions of equilibrium for 
static and rotating stars and some phenomenological assumptions. The resulted approximate 
mass and radius formulae are compared with data obtained from numerical integration of 
gravitational field equations in case of neutron stars described by realistic equations of state 
and strange stars described by the bag model equation of state. 

The present paper is organized as follows. The general formalism allowing to obtain mass 
and radius formulae for rotating general relativistic stars is presented in Section 2. In Section 
3 we apply our results to the case of neutron stars described by realistic equations of state. 
In Section 4 we consider the case of strange stars. We discuss and conclude our results in 
Section 5. 



II. THE GENERAL FORMALISM 

For a static equilibrium stellar type configuration, with interior described by the metric 

ds 2 = e 2 ^dt 2 - e x ^dr 2 - r 2 (d6 2 + sin 2 6d<p 2 ) , (5) 

the condition of the hydrostatic equilibrium, which follows from the Bianchi identities, can 
be written as 

2z/ + 2 / = const., (6) 

J P + P 

with p and p the energy density and the pressure of the matter respectively ( in the present 
paper we shall use units so that c — G—1). 

We shall also make the assumption that the matter EOS is a one parameter dependent 
function, 

p = p(n),p = p(n), (7) 

with n the proper baryon density. For such an equation of state the heat function defined 
by 

H{n)= i 1 m, (8) 

v ' J p(n) + p(n) dn w 
is a regular function of n.It can always be written in the form [16] 

H(n) = hi f(n), (9) 

where 

/(») = ^1±M, (10) 

is the enthalpy per baryon. 

As applied at the center and at the surface of the star respectively, the hydrostatic 
equilibrium condition yields 
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(2i/ + 21n/(ra)) | center = (2z/ + 2 In f(n)) \ sur f ace . (11) 

At the vacuum boundary of the static star the Schwarzschild exterior solution gives the 
metric and we have [17] 

2^ = ln(l-^V (12) 



Rstat 

with M sta t and Rstat the mass and radius of the static stellar configuration. We denote by 
nc and n$ the baryon density at the center of the star and at the surface, respectively. 
Therefore from ( |TTD and flT2] ) we obtain the following general and exact expression for the 
mass-radius ratio of the static star: 

M stat 1 



Rsinf. 2 



[l-C{nc)F(nc,n s )], (13) 



where we have denoted 



and 



F(n c ,n s ) = p^, (14) 



C (n c ) = exp (p c ) = 9oo \r=o ■ (15) 



For a given equation of state C (nc) is a function of the central density only. From a 
physical point of view C (n c ) can be related via the relation ~ — 1 = z c to the redshift z c of 
a photon emitted from the center of the star. 

For the mass of the star we can obtain another general representation by assuming 

M stat = ^^RLt A (n c ) , (16) 

where p spec is a specific density that can be arbitrarily chosen (for example it is the central 
density of the minimum mass configuration) and A (nc) is a function describing the effects 
of the variation of the central density on the basic parameters of the stellar configuration. 

The two unknown functions C (nc) and A (nc) can be determined by fitting equations 
(O) and (ITB) with the exact values of the mass and radius obtained by numerically in- 



tegrating the TOV equation for a given equation of state. Their knowledge allows us to 
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construct exact mass and radius formulae for sequences of static general relativistic stars 
having different central densities. From equations fll^D and fllBj ) we obtain the radius and 
the mass of the star in the form 

/ o i \i/2 

^=(lSI<^ \l-Cinc)F(n c , ns) r, (IT) 



4vr . , , / 3 



\ 3/2 

M stat = ^ Pspec A (nc) ( ^r-rr-^ t 1 ~ c ^ F ^ ns ^' 2 ■ ( 18 ) 

3 \167T A(n c )p spec J 

Solving the equations d ^^ Qi = and = for the value of the central density n c 

will give, with the use of fllTD and (p~8[ ) the maximum values of the radius and mass of the 
static star, respectively. 

To describe the interior of the rotating general relativistic star we shall adopt the for- 
malism presented in [16] and [9]. Under the hypothesis of stationarity, axial symmetry and 
purely azimuthal motion a coordinate system (t, r, 9, 0) can be chosen so that inside the star 
the line element takes the form 

rfs 2 = ^2^2 _ B 2 r 2 gin 2 Q ^ _ _ A 2 ^2 + ? (1Q) 

where N,N^,A and B are functions of r and 9 only. As measured by the locally non- rotating 
observer the fluid 3- velocity is given by U — Br sin 9^^- , where = % is the angular 
velocity of a fluid element moving in the ^-direction (physically it is the angular velocity as 
measured by an observer at spatial infinity) [9], [16]. 

From the point of view of our phenomenological approach the most important result is 
the equation of the stationary motion, which results from the Bianchi identities and which 
for a rotating perfect fluid reduces to [16] 

1 *+£-A(i»r) = -s?l (20) 



p + p dx l dx i dx % dx l 

where S = T 2 ^Ursin9 , v = InN and T = (l~U 2 y 1/2 . If Q = const, (case called 
uniform or rigid rotation) equation ( P0| ) can be integrated to give the following fundamental 
result describing the stationary equilibrium of a rotating general relativistic star: 



(2v + 2\nf(n) +ln(l - U' x 



I center 



(2u + 2\nf(n) +ln(l - f/ : 



I surface 



(21) 



Equation fl2"T|) is just the generalization to the case of rotation of the well-known static 
Bianchi identity 2u + 2 J = const. [17] we have already used to describe static stellar 
configurations. 

We assume that the vacuum boundary of the rotating star is described by the Kerr metric 
[1] in the Boyer-Lindquist coordinates, which has the form [18]: 

2M rot r \ ( r 2 + a 2 cos 2 9 



ds' 



r 2 + a 2 cos 2 9 



dt 



(r 2 + a 2 sin 2 6)d6 2 - 



\r 2 + a 2 -2GM rot r 
r 2 + a 2 ) sin 2 9 + 



r 2 + a 2 cos 2 9 dr z 



(22) 



2M rot ra 2 sin 4 9 



r 2 + a 2 cos 2 9 



2 2M ro jrasin 2 # 

b ~ 2 2 i 2 2a dtd 4- 

r z + a z cos z t7 



In this form the Kerr metric is manifestly axially symmetric and closely resembles the 
Schwarzschild solution in its standard form.M rot is the mass of the source and the parameter 
' is the ratio between the angular momentum J and the massM rot of the rotating 



Mrot 



star. 



Let's apply equation (|2~ID at two points: at the center of the dense core and at the pole 
of the rotating star. We denote by C^t the value of the metric tensor component g 00 at the 
center of the star. At the polar point 9 = and r = R p = const., where R p is the polar 
radius of the star. We also have U = . Therefore at this point the line element is given by 

2M rot 



ds loi 



R„ 



(l + a 2 /R 2 p )~ l 



dt 1 



(23) 



Consequently from equations ( pT|) and ( p3|) we obtain the following exact mass-polar 
radius relation for the rotating general relativistic configuration: 



M, 



rot 



R p 2 



(l + a 2 /R 2 p ) 1 [l - C« (no, •••) F (n c , n s ) 



(24) 



The function F (nc, ns) is identical to that of the static case. Let's apply now equation 
(PT|) for two points situated in the equatorial plan of the rotating star: at the center of 
the star and at the equator respectively. At the equator r = R e = const, and 9 = | 
(R e is the equatorial radius of the star). For a uniform rotation the rotation angle of the 
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source/observer at the equator is = Qt . Taking into account these results we obtain the 
Kerr metric at the equator of the rotating star in the form: 



2M 



rot 



dt 2 . 



(25) 



With the use of fl2"ip and (E5I) we obtain the following mass-equatorial radius relation: 



M rot = 1 

Re 2 

where a new function 



1 - eg (no, Q, R e ) F (n , n s ) - R 2 e Q 2 (l + ^) 



(26) 



C% (no, 0, i2e) = (no, 0, i? e ) / ( 1 - U 2 ) \ r=R 



(27) 



has also been defined. 



From Eqs. (53) and ( p6|) we obtain the ratio e of the polar and equatorial radius (the 
oblateness) of the star in the form: 

e= Rp = 1 + jj 1 ~ F ("C "g) ~ ^ j 1 + j) (2g) 

e ~Re~(l + aQy l-C^(n c ,n,R e )F(n c ,n s ) 
In the case of the static star we have also proposed the alternative Eq.([TJ]) for providing 
another mass-radius relation. We shall generalize this equation to the rotating case by 
assuming that the following formula relates the mass of the rotating star to its equatorial 
radius: 



M ra t — ^ P RlCrrt (no, 0, Re) 



(29) 



with C^ol (no, ^, -R e ) a function describing the combined general relativistic effects of rotation 
and central density on the mass of the star and with p spec the same specific density as used 
in the static case. For Q = 0, R = R e and Eq.(|29|) must reduce to the static case equation 

Equations (pop , (|28|) and ( |29| ) give a complete and exact description of the mass and 
radius of the rotating general relativistic star. But unfortunately the present approach, 
which is basically thermodynamic in its essence, can not predict the exact form and the val- 
ues of the three unknown functions entering in the formalism. The only thing we can 



do is to assume, also based on the static case, some empirical forms for the functions 
Crot (nc, R e ) ,i = 1,2,3 and to check whether the resulting formulae can give a satis- 
factory description of rotating star configurations. Therefore in the following we shall use 
the following five approximations: 

i) As a first approximation we shall define the moment of inertia I rot of the rotating 
compact relativistic object via the Newtonian expression 

I rot = \M rot R\. (30) 

In fact over a wide range of M rot and R e the corrections added by general relativistic effects 
to the moment of inertia can be approximated by I rot = 0.21 Mr 2 °M^ e t [15], but we shall not 
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use this representation. By adopting the Newtonian formula we obtain 



tt =-£=-^ fi ' < 3i) 

ii) We assume that the function C^t ( n c -Re), describing the metric tensor component 
e u at the center of the rotating star is given by 

cffl (n c , n, R e ) = [(i-\g (^) c (n c ) , (32) 

with g (^77) an EOS dependent function given by 

(M stat \ a 

tlstat 

and a a non-negative constant. C (n c ) is also the function corresponding to the static case. 

iii) We assume that inside the rotating star C^t (nc, Re) is independent of the an- 
gular velocity Q of the rotating compact object and can be represented by the function 
corresponding to the static case: 

cl 2 J t (n c ,{l,Re)=C(n c ). (34) 



iv) We suppose that 



Clf t (n c ,n,R e )=A(n c )e. (35) 
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A (tic) is again the function corresponding to the static case. 

With these four phenomenological assumptions Eqs. (|26|)-(|29|) lead to the following 
representation of the basic physical parameters of the rotating general relativistic star: 

-1/2 



D / stdt \ 



R 



stat 



M stat a^ 2 



\ L + 25 e2 


) 


- i^ 2 (l + !^ 2 ) 




(l - f Bff!^ 


2 

1- 


< j 2M stol > 




1 2 1 iM stat ) 



(36) 



(37) 



M rot = ^P^lA (n c ) Rle = M stat ( £-) * e. (38) 



3 V -Rstat / 

We shall also suppose that a is a universal constants and we shall choose a ~ 5. In this 
formulation of the general relativistic problem of the rotation the oblateness parameter e of 
the star is given by the roots of the third order algebraic equation fl37|) . 

The equatorial radius is defined only for values of the angular velocity satisfying the 
condition ~ Q 2 < R ^f at .Therefore for the maximum admissible constant angular velocity of 

4 H stat 

the maximally rotating star in uniform rotation we obtain the relation 

A/rstat \ / Tjstat \ 

max 1 / ^"max 1 

M ) \10km ) 



V<C , (39) 



2 (M l 1 ^ 2 

C = ^VG K ®; = 10330.489s- 1 , (40) 



where 



/n 10 9 

for a = 5. 

Equation ( |3~9"D is very similar to the "empirical" formula discussed in [l]-[7]. The co- 
efficient of proportionality in (^) is independent of the equation of state of the dense 
matter, but its numerical value does not fit the calculated value. On the other hand for 

/ Turstat \ 1/2 / nstat \ —3/2 

Q = C [^j^J (lofc^J the radius of the star tends to infinity. 

An alternative expression can be obtained by imposing the restriction 
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/ TUT rot ICtr M stat I M stat 

pmaxQ a / lvl max ^ I ^Al lvl m g X ^ / lvl max I a -i \ 

n e "max ~ P ~ \l jj rot ~ A/ r R stat ~ \/ nested ' ( " ' 

L max V L max V max 



where /3 is the maxium allowable equatorial speed of the star and we have also used the 
Newtonian force balance equation between the gravitational and centrifugal force. Therefore 
we obtain 



V 2 / ostat \ - 3 / 2 



o ~ n 

J L max ~ " 



max 1 / max 1 

^ M / / 



4 n#sta< 

_|_ m<i:r 



2 f?stai 

max 



-1/2 



(42) 



It is interesting to note that the values of 2M^ f x / R^ x are in a narrow range of (0.467, 
0.667) [8], [12]. Taking for 2M^/Rg*. a mean value of 0.58 [12] it follows that (3 « 0.54 
and we find 

tU.„ 0s ^j (iiS) ■ (43) 

with Cs = 7708s -1 , value which coincides with that proposed in [8] and differs only within 
1% from the value Cs = 7840s -1 obtained in [9]. 

By taking for twice of the ratio of the maximum static mass and radius a mean value of 
0.58 [12] and considering (3 ~ 0.54, we obtain 

/ Afstat \ V 2 / r>stat \ ^ 3 / 2 

a^«Cs[-^) (^j , (44) 

with Cs = 7708s -1 , value wich coincides with the value obtained in that differs within 5% 
from the value Cs = 7840s -1 obtained in [20]. 

The maximum radius of the maximally rotating configuration can be obtained from 



R ™™ ~ ?W and is given by 



l Dstat 
\ *"max 

where 



/ Ajstat \ 

ryrot ^ si I ±y± max \ jDStat (AK\ 
''max 



-1/2 



lOfcmCsVlO^mV^; ' 1 ] 

The mass M™* xrad of the maximal radius rotating neutron star follows from fl58p and is 
given by 
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' yrstat \ 

nfrot = C l max 1 jystat ( A7\ 

maxrad ^ RM \ jw £ I lvl maxi \ ) 



where 



R s 

- 1 v max 



/ ±v± max I /^y3 | max I ^rac( | max \ ( AQ\ 

~DstaT J ~ C R\ ~DstaT 6 max "Tw^T > l 48 J 
\ max / \ max / \ max / 



e 



rad I max 

™™ I JDstat 

\ L max , 



1 


i-(i-/5 2 ) 7 (1- 


2M^«n 

f?stat } 
± *"m ax / 


/? 2 (1 + ^ 





(i + r; 


2 i (i 2 ^r) 




\ ^^max / _ 





(49) 



The maximum mass M™* x °f the rotating star can be obtained from the equation 



dM M stat d , 3 x M stat /„ p2 ai2e R3 <9e\ 

of |n=fW ~ V^ e J |n = n — - ^ e + ^Mj |n=n — - °' (50) 

leading to 

ro ,_ c(^,iu) _ ^te)^ 

i,i„,x { il max) D3(n \ #3 (O \ ' ^ ' 



where Cm ( ^f 3tat ) is a dimensionless EOS dependent function. Therefore for the maximum 
mass of the maximally rotating configuration we obtain 



MZL = C U (^)M«£. (52) 

V -ttstat / 

Equations ( fiEf) and (p2|) show the existence of a proportionality relation between max- 
imum mass and radius of the rotating and non-rotating configuration, respectively, as has 
already been suggested, on an empirical basis, in [12]. 

An investigation of 12 EOS performed in [12] for 12 realistic EOS of nuclear matter led 
to a (mean) value of Cm = 1.1807 , while the same calculation performed by us with the use 
of data presented in [20] for other 14 different EOSs gives Cm = 1-177 , leading to a mean 
value of Cm = 1-1790 for the 26 considered EOSs. 

Eqs. ©-© lead, for (3 w 0.54 and 2M^fjR^ x w 0.58 to values of e max = 0.365 and 
Cm = 1.233. 

Therefore we may conclude that the proportionality between the maximum mass of the 
rotating star and the maximum mass of the static configuration is universal, being with 
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a very good approximation independent of the equation of state of dense matter. On the 
other hand the relation between the radius of the rotating and non-rotating configuration 
is EOS dependent, the coefficient of proportionality slightly decreasing with the increase of 
the mass-radius ratio of the static star. 



III. APPLICATIONS TO NEUTRON STARS 

The correct mathematical and physical modelling of millisecond pulsars can be done 
only in the framework of general relativistic equilibrium models for rapidly rotating neutron 
stars. Such models are solutions of the Einstein's equations for the axisymmetric stationary 
gravitational field and they must be constructed numerically. Recently several independent 
numerical codes have been developed by different groups of researchers and have been used 
to obtain rapidly rotating neutron star models based on a variety of realistic equations of 
state. Hence a large amount of numerical data is now available. 

In the present Section we shall apply the results of the phenomenological formalism 
presented in the previous Section to the case of neutron stars described by realistic equations 
of state. The data are selected from the paper by Cook, Shapiro and Teukolsky [20], who 
constructed general relativistic rotating star sequences for 14 nuclear matter equations of 
state. Detailed data are presented only for 5 equations of state. For the sake of comparison 
we have used the equations of state denoted A [21], AU [22], FPS [23] and L [24]. At low 
densities all these equations of state employ the Feynman, Metropolis and Teller [25] EOS 
and then join onto the Baym,Bethe and Pethick [26] EOS up to neutron drip. The equations 
of state are given in a tabular form and small changes in the way the tabulated equation of 
state is constructed do have a small effect on the resulting neutron star model. 

In Tables 1-4 we present the comparison of the basic physical parameters of rotating 
stars obtained, for these four equations of state, with the help of Eqs. (|36"|)-(|3"8|) and by 
numerically integrating the gravitational field equations. 
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TABLES 



nis- 1 ) 


Rnum{km) 


Rcalc(km) 


Rnum 


M num [Mo 


M calc (M & ) 


M ca lc 





9.586 


9.586 


1.0 


1.40 


1.40 


1.0 


3244.1 


9.763 


9.899 


1.014 


1.4030 


1.423 


1.014 


5018.9 


10.06 


10.393 


1.033 


1.4077 


1.456 


1.034 


6136.6 


10.38 


10.876 


1.047 


1.4123 


1.483 


1.050 


6940.0 


10.74 


11.340 


1.055 


1.4169 


1.506 


1.063 


7544.8 


11.14 


11.777 


1.057 


1.4214 


1.525 


1.073 


7953.7 


11.56 


12.125 


1.048 


1.4252 


1.538 


1.079 


8236.4 


12.00 


12.397 


1.033 


1.4285 


1.547 


1.083 


8431.1 


12.49 


12.602 


1.008 


1.4313 


1.554 


1.085 


8590.6 


13.72 


12.781 


0.931 


1.4340 


1.559 


1.087 



TABLE I. Comparison of the radius R ca ic an d mass M ca i c obtained from Eqs.(36)-(38) with the 
radius R n um and mass M num obtained by numerically integrating the gravitational field equations 
[20] for EOS A [21]: 1.4Mq normal mass sequence 
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/?,,„/^f km) 




- 1 VJ - nuin \ - 1 V - L y ) 


Mr-nlAMr^ 

xvx ca,lc\- LV - L (i) } 


M ca .lc 





10.85 


10.85 


1.0 


1.4 


1.4 


1.0 


2883 8 


11.08 


11.259 


1.016 


1.4031 


1.424 


1.015 


4112.1 


11.36 


11.736 


1.033 


1.4065 


1.449 


1.030 


5033.6 


11.71 


12.268 


1.047 


1.4103 


1.472 


1.044 


5712.1 


12.10 


12.791 


1.057 


1.4140 


1.491 


1.055 


6174.6 


12.49 


13.233 


1.059 


1.4173 


1.504 


1.061 


6544.0 


12.94 


13.649 


1.054 


1.4206 


1.515 


1.066 


6841.8 


13.51 


14.035 


1.039 


1.4239 


1.522 


1.069 


7016.9 


14.06 


14.286 


1.016 


1.4263 


1.527 


1.070 


7165.0 


15.45 


14.519 


0.939 


1.4287 


1.530 


1.071 



TABLE II. Comparison of the radius R ca ic an d mass M ca i c obtained from Eqs.(pq)-(|3q) with the 
radius R n um and mass M num obtained by numerically integrating the gravitational field equations 
[20] for EOS FPS [22]: 1.4M normal mass sequence 
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( krri} 


Rnurn 


A v± niiTYi V W / 


- Ly - L CQ,lC\- LV - L (i) ) 


M ca .lc 





9.411 


9.411 


1.0 


2.1335 


2.1335 


1.0 


3784.6 


9.744 


9.669 


0.992 


2.1417 


2.174 


1.015 


5827.1 


10.05 


10.060 


1.001 


2.1547 


2.216 


1.028 


7381.3 


10.40 


10.525 


1.012 


2.1703 


2.304 


1.061 


8507.7 


10.78 


10.987 


1.019 


2.1863 


2.371 


1.084 


9308 9 


11.20 


11.403 


1.018 


2.2018 


2.430 


1.104 


9854.5 


11.64 


11.740 


1.008 


2.2157 


2.477 


1.118 


10211.0 


12.10 


11.988 


0.990 


2.2278 


2.512 


1.127 


10426.0 


12.57 


12.150 


0.966 


2.2373 


2.534 


1.132 


10587.0 


13.66 


12.279 


0.898 


2.2467 


2.551 


1.135 



TABLE III. Comparison of the radius R ca ic an d mass M ca i c obtained from Eqs.(|36|)-(|38|) with 
the radius R nU m and mass M num obtained by numerically integrating the gravitational field equa- 
tions [20] for EOS AU [23]: maximum mass normal mass sequence 
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( hrn] 




- 1 VJ - num \ - 1 V - L y ) 


- LV± CCLtC\ ±yj -(i) } 


M ca .lc 





13.70 


13.70 


1.0 


2.7002 


2.7002 


1.0 


2212.2 


14.20 


14.01 


0.986 


2.7063 


2.739 


1.012 


3495.1 


14.68 


14.52 


0.989 


2.7181 


2.800 


1.030 


4485.0 


15.24 


15.13 


0.993 


2.7331 


2.869 


1.050 


5212.1 


15.88 


15.75 


0.992 


2.7494 


2.936 


1.068 


5712.7 


16.57 


16.29 


0.983 


2.7650 


2.991 


1.081 


6051.2 


17.30 


16.72 


0.966 


2.7797 


3.033 


1.091 


6263.5 


18.05 


17.02 


0.942 


2.7923 


3.062 


1.096 


6422.3 


19.14 


17.26 


0.901 


2.8056 


3.086 


1.099 


6482.9 


20.66 


17.36 


0.840 


2.8126 


3.095 


1.100 



TABLE IV. Comparison of the radius R ca ic an< l mass M ca \ c obtained from Eqs.(|36|)-(|38|) with 
the radius R nU m and mass M num obtained by numerically integrating the gravitational field equa- 
tions [20] for EOS L [24]: maximum mass normal mass sequence 
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For 1.4M sequences described by EOS A and FPS the maximum error of our prediction 
is around 9%. For the maximum mass normal sequences of EOS AU and L the maximum 
error in the predicted value of the mass does not exceed 14% but for large angular speeds it 
is around 16% for the radius of the rotating compact object. 

IV. APPLICATIONS TO STRANGE STARS 

It is generally believed today that strange quark matter, consisting of u-,d- and s quarks 
is energetically the most favorable state of quark matter. Witten [27] suggested that there 
are two ways of formation of the strange matter: the quark-hadron phase transition in the 
early universe and conversion of neutron stars into strange ones at ultrahigh densities. In the 
theories of strong interactions quark bag models suppose that the breaking of physical vac- 
uum takes place inside hadrons. As a result the vacuum energy densities inside and outside 
a hadron become essentially different and the vacuum pressure on a bag wall equilibrates 
the pressure of quarks thus stabilizing the system. 

If the hypothesis of the quark matter is true, then some of neutrons stars could actually 
be strange stars, built entirely of strange matter [28], [29]. Caldwell and Friedman [30] have 
presented arguments against the existence of strange stars. For a recent review of strange 
star properties see [31]. 

There are several proposed mechanisms for the formation of quark stars after galaxy 
formation. Strange stars are expected to form during the collapse of the core of a massive 
star after the supernova explosion as a result of a first or second order phase transition, 
resulting in deconfined quark matter [32]. Another possibility for strange star formation is 
that some rapidly spinning neutron stars in low-mass X-ray binaries (LXMBs) can accrete 
sufficient mass to undergo a phase transition to become strange stars [33]. In this scenario 
it is supposed that at the beginning of accretion the mass of the neutron star is 1.4. It has 
been shown [34] that the amounts of matter accreted by 18 millisecond pulsars in binary 
systems exceed O.5M . Hence some of the millisecond pulsars may be strange stars. Strange 
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stars have also been proposed as sources of unusual astrophysical phenomena, e.g. soft 7-ray 
repeaters [35] , pulsating X-ray burster [36], cosmological 7-ray bursts [35], [37-38] etc. The 
mechanism of the phase transition from neutron to quark stars in low LXMBs also results 
in the excitation of stellar radial oscillations that can be damped by gravitational wave 
radiation instead of internal viscosity [39] .The discovery of kHz quasi-periodic oscillation in 
LXMBs [40] implies that the compact stellar object must have very soft equation of state, 
which is consistent with that of strange stars [41-42]. 

Assuming that interactions of quarks and gluons are sufficiently small the energy density 
p and pressure p of a quark-gluon plasma at temperature T and chemical potential \if (the 
subscript f denotes the various quark flavors u,d,s etc.) can be calculated by thermal theory. 
Neglecting quark masses in first order perturbation theory and supposing that quarks are 
confined to the bag volume (in the case of a bare strange star, the boundary of the bag 
coincides with stellar surface), the equation of state is 

P ^^, (53) 

where B is the difference between the energy density of the perturbative and non- 
perturbative QCD vacuum (the bag constant). Equation fl53] ) is essentially the equation 
of state of a gas of massive particles with corrections due to the QCD trace anomaly and 
perturbative interactions. These are always negative , reducing the energy density at given 
temperature by about a factor two [43]. In the limit p — > (at the star's surface) we have 
p — > AB. The equation of state (|53|) does not depend upon quark flavor number, hence it will 
be correct either for strange quark matter (m s — > 0) or for normal quark matter (m s — > 00). 
For any intermediate values of m s the state equation (^) gives the pressure with error 
less than 4% [31]. Thus the equation of state of strange matter is mainly determined by the 
vacuum energy density B. 

The bag model equation of state ( |53"D has been the basis for the study of most of the static 
relativistic models of strange stars [27], [29]. Based on the numerical integration of the mass 
continuity and hydrostatic equilibrium TOV (Tolman-Oppenheimer-Volkoff) equations for 
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different values of the bag constant these authors obtained a complete description of static 
strange stars. Using numerical methods Witten [27] and Haensel et al. [29] obtained the 
maximum gravitational mass M max , the maximum baryon mass Ms,max = 1-66 x I0~ 27 kg x 
Ng (iV^-the total baryon number of the stellar configuration) and the maximum radius R m ax 
of the strange star ,as a function of the bag constant, in the form [27], [29], [9]: 

1.9638M© 2.6252M 10.172A;m 

^max I p j max 1 p > ^max 1 p j v^^v 

V -D60 V -D60 V -D60 

where B 60 = B / \60MeV fm~ 3 ). 

Colpi and Miller [44] and Glendenning and Weber [45] have investigated the rotational 
properties of strange stars in the slow rotation approximation. As far as rotational defor- 
mations are concerned, there are a number of detailed differences between the strange star 
models and standard neutron stars. Exact numerical calculations of rapidly rotating strange 
stars were done by Lattimer et al. [10], Gourgoulhon et al. [9] (by using a multi-domain spec- 
tral method that enable to treat exactly the density discontinuity at the surface of strange 
stars) and by Stergioulas,Kluzniak and Bulik [42]. Rotation increases maximum allowable 
mass of strange stars and the equatorial radius of the maximum mass configuration. Gour- 
goulhon et al. [9] obtained for the maximum mass and radius of quark stars the following 
two exact formulae 



_ 2.831M 16.54fcm 

m max - fn— ' U eq,M max ~ 7W= • I 00 ) 

V-D60 V-D60 



In the present section we shall derive, by using the formalism presented in Section I, 
analytic mass and radius formulae for general relativistic static and rotating equilibrium 
strange matter configurations described by the bag model equation of state (|53|).We shall 
begin with the study of the static strange star, but presenting and alternative and physically 
more involved discussion of this case. 

The changes caused by the general theory of relativity in the conditions of thermal 
equilibrium,taking into account the gravitational field of the bodyare of fundamental im- 
portance. In a constant gravitational field we must distinguish the conserved energy Eq of 
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any small part of the stellar object from the energy E measured by an observer situated 
at a given point. These two quantities are related by E = E^fg^ [46], where goo is the 
time component of the metric tensor. A similar change occurs in the condition of the con- 
stancy of the chemical potential throughout the star. The chemical potential is defined as 
the derivative of the energy with respect to the number of particles N, p — (^)s,v- Since 
this number is a constant for the stellar object, N = constant, for the chemical potential 
measured at any point inside the gravitating body we have the relation [46]: 

fiyfgoo = constant. (56) 



A similar relation also holds for the temperature T, T^fg^ = constant., since we suppose 
that the strange star is in thermal equilibrium [46]. Consequently, ^ = constant inside 
the compact object. Hence = ^p-.At constant volume (equal to unity) we have dp = 
sdT + ndp,where s and n are the entropy and number of particles in unit volume of the 
bodyrespectively. With the use of dT = T— and taking into account that pn + sT = p + p 
we obtain the following equation relating the equilibrium chemical potential to the energy 
density and pressure of the star [46]: 

dp = dp 
p p + p 

Consider now a static equilibrium quark matter configuration satisfying the bag model 
equation of state fl5"3|) . Let us compare the values of the chemical potential p = A 4 / a t two 
points:at the center of the star and at the vacuum boundary. From equation (|56|) we obtain 

[^V^OO] center = IP'S V^Oo] boundarv , (58) 



J center " i»UUJ boundary 

where the indices C and S refer to the center and to surface of quark star respectively.At the 
vacuum boundary the gravitational field of the strange star is described by the Schwarzschild 
solution, which gives [17]: 

2M stat 

9oo = 1 - ~b , (59) 

r£.„ 



i-stat 
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where M sta t and R s t a t are the total mass and radius of the static strange star,respectively.At 
the center of the star the time component of the metric tensor has a constant value (this 
also follows from the Bianchi identity J:ln(<7oo) = — [17] an< ^ we denote 

#00 \center= C (pc , B) . (60) 

From a physical point of view C(pc, B) can be related via the relation ^ — 1 = z c to the 
redshift z c of a photon emitted from the center of the quark star. For a given static strange 
matter configuration the value of C depends only on the central density of the quark star 
pc and on the bag constant. Therefore from equation fl58f ) we obtain 



&= 1 fi-gj***) (61) 
4 C{ PC ,B) \ R stat > 



With the use of the bag model equation of state ([53]) we can integrate equation (|57j) to 
obtain 

/i = C (p- J B) 1/4 . (62) 

The integration constant Cq can be determined by calculating the chemical potential at 
the center of the quark star. Hence we obtain 

The variation of the chemical potential inside the quark star can be represented as 



At the surface of the star p$ — > AB. Therefore from equation fl64j) it immediately follows 
that 

-- (3B)1/4 (^V- (65 » 

In order to simplify notation we shall introduce a dimensionless parameter rj = 
so that C(pc,B) = C(rj) .By eliminating ^ from equations (pi]) and (|5"5|) we obtain the 
following exact formula for the mass-radius ratio of a strange star: 
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M. 



stat 



(66) 



Rstat 

For a given equation of state the mass-radius ratio of the star depends on the values of 
the metric tensor component at the center of the star, £7(7?), only. A possible representation 
for the function giving the values of #oo at the center of the quark star is in the form of a 
power series C(rj) = const. + ^with a>i constants. 

As applied on the star surface the mass continuity equation leads to a rough approx- 
imation of the quark star mass of the form dM ^ at ~ ^ stat ~ lQnBR 2 stat . A mass-radius 
relation of this form could also describe zero pressure quark matter, with e = AB,p = 
and M s tat(Rstat) = W7TB ^ stat . But for densities greater than AB the effects determined by 
the large central density become important. Hence for strange quark stars we propose the 
following mass-radius relation: 

Mstat = -^-a( V )R 3 stat , (67) 

with 0(77) a function describing the variation in the quark star mass due to the increase of 
the central density. 

The exact form and the values of the functions C(rj) and a(n) can be determined only 
by numerically integrating the gravitational field equations. By fitting the numerical data 
given in [47] for the mass and radius of the strange star with the expressions (|66|) and (|67|) 
we obtain the following representations for these functions (in the present paper we shall 
consider B = 10 1A g/cm 3 = 56Mev/ fm 3 ): 

C(rj) = 44.005^- - 6.68158^r + 2.7403- + 0.0554667, (68) 
a{rj) = 0.0000521833?7 3 - 0.00378523// 2 + 0.114564r/ + 0.624094. (69) 



The numerical constants in equations (|6lf ) and ( |69"D depend on B because the numerical 
data have been calculated at a given B. For the polynomial fittings (^) and ( pP| ) the 
correlation coefficient r = 0.9997 and the probability P < 0.001. Therefore for a given value 
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of the bag constant B we obtain the following exact representations for the radius and mass 
of the static strange matter configuration obeying the MIT bag model equation of state: 

RmM = (_i£_) 1/2 (_L L _ c M {^f))\ (to) 



The variations of the radius and mass for a strange star (B = lO 1 ^ /cm 3 ) as a function of 
the parameter rj are represented in Figures 1 and 2. For the sake of comparison we have also 
presented the data obtained by numerically integrating the TOV and hydrostatic equilibrium 
equations [4 7]. Using (|6"8"D-(f7lD we can reproduce the values of the mass and radius of the 
quark star obtained by numerical integration with an error smaller than 1%. The maximum 
radius of the strange star is obtained from the condition ^ = O.The corresponding algebraic 
equation has the solution r)^ ax = 9.99012 (this value depends of course on the value of 
B), giving the value of the ratio of the central pressure and bag constant for the maximum 
allowable radius R max of the static strange star. This can be expressed as 



R%L = °- 569906 x V^g' ( 72 ) 

and its numerical value for B = 10 1A g/cm 3 is R m ax — 1-1436 x 10 6 cm.From the condition 
^ = it follows that r/j^l = 22.41173 and the maximum mass of the static quark star is 
given by 



3 

r .,/,,/ ™r,o^ . V 3 / C 



JIC = 0.297866 x r - . (73) 

2(32tt5)5 \GJ k j 

From equation (|73| ) and for the chosen value of the bag constant, we obtain a value of 
M max = 2.016M Q .These results are in good agreement with the previously proposed (Witten 
[27],Haensel,Zdunik and Schaeffer [29]) maximum radius and mass values,given by equations 
(p4j) (from equations ( [54] ) and for B = 56MeV/ fm 3 we obtain M max = 2.O3M ).For values 
of i] > static quark star models would be unstable to radial perturbations. 
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As an application of the mass and radius formulae obtained for the static strange stars 
we shall derive an explicit expression for the total energy of the quark star. The total 
energy (including the gravitational field contribution) inside an equipotential surface S can 
be defined, according to Lynden-Bell and Katz [48] and Gron and Johannesen [49] to be 

E = E M + E F = / [K]ds, (74) 

<57T JS 

where is a Killing vector field of time translation, £ s its value at 5* and [K] is the jump 
across the shell of the trace of the extrinsic curvature of S, considered as embedded in 2- 
space t = constant. Em = SsTi£, l \/—gdSk and Ep are the energy of the matter and of the 
gravitational field, respectively. This definition is manifestly coordinate invariant. In the 
case of the static strange star with the use of equation (|70|)and (fTl|) we obtain for the total 
energy (also including the gravitational contribution) the following exact expression: 

1/2 

(75) 



E = E SQM + E F = -J -^Cirj) 



32ttB 



1 "- 1 / 1 _ C() , ) ,M 



a(n) 3 y V 3 

where E S qm is the total energy of the quark matter. 

The variation of the total energy of the strange star as a function of the parameter n is 
represented in Figure 3. 

The minimum value of the total-matter plus gravitational-energy of the strange matter 
configuration is obtained for f] m i n = 5. 68171. The most stable static stellar configuration 
made of strange matter is given by quark stars with radius R s tab = 9.97179 x 10 5 cm and 
with mass M stab = O.96558M , corresponding to values of the central density of the order 
of p c = 5.681 x B = 5.681 x 10 u g/cm 3 . 

We shall consider now the study of the rotating strange star configurations. We shall 



compare our results obtained with the use of equations (|36|)-(|38|) and (|70D-(|71~D with the 
results provided by Stergioulas,Kluzniak and Bulik [42] and obtained by numerically in- 
tegrating the gravitational field equations for maximally rotating ("Keplerien") models of 
strange stars. The results of Stergioulas et al.[42] are also in very good agreement with the 
results of the exact numerical models of rotating strange stars built of self bounded quark 
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matter of Gourgoulhon et al. [9], the difference between these two works being smaller than 
1%. In order to improve the accuracy of the expressions (^)-(^) we shall consider that the 
function C^ t (77, R e , can be expressed in a more general form as 

C^ t ( V} R e ,n) = (l-R 2 e n 2 Yc( V ), (76) 



and we will assume that the parameters a and 7 are not constants, but some angular 

velocity dependent functions given by 

1.23188 x 10 9 2.15445 x 10 5 
a = — + 13.1399, (77) 

7 = 9.66251 x 10 U Q- 3 - 4.56952 x 10 8 1T 2 + 6.99642 x lO 4 ^ 1 - 2.54401. (78) 



Equations (|77|)-(|78|) take into account the variation of the central density of the maxi- 
mally rotating strange star due to the increase of the angular velocity. 

In Fig. 4 we have represented the variation of the radius of the strange stars, given by 
equation (|79|) together with the angular velocity dependent functions a and 7 and the 
values given in Stergioulas et al.[42], calculated for the same values of the central density 
and angular velocity.The mean of the difference between these two sets of values is smaller 
than 1%. The variation of the radius of the maximally rotating strange star as a function 
of both central density and angular velocity is represented in Fig.5.Figs.6 and 7 present the 
variation of the mass of the rapidly rotating strange star as a function of the angular 
velocity Q and of Q and central density, respectively. 

From equation (^) and equations (|70|) and ([H]) it follows that the radius of the rotating 
strange star can be expressed, as a function of the central density and angular velocity 
onlyin the following form: 



Re 



12 i-cWV 



(79) 



\ 647rBa(rj) - 3afl 2 ' 
Hence in this approximation the basic rotational parameters of maximally rotating 
strange star can be represented in terms of the static configuration and of the angular 
velocity only. 
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V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have suggested the possibility of the existence of a universal 
pattern expressing the basic properties of rotating compact object as simple functions of the 
parameters of the static object and of the angular velocity only. We have obtained exact 
formulae which give the dependence of the radius and mass of the static and rotating stars 
on the central density of the stellar object and of its angular velocity.In the static case this is 
made possible due to the constancy of the chemical potential. The two unknown functions 
involved in the model must be obtained by fitting the exact formulae with data obtained 
from the numerical integration of the structure equations of the neutron or quark star. The 
resulting analytical expressions can reproduce the radius and mass of the strange star with 
an error smaller than 1% and they also provide a simple way to obtain the maximum mass 
and radius of the static configuration. 

In the rotating case, with the use of the hydrostatic equilibrium condition, which is 
the consequence of the Bianchi identities, we have also obtained exact mass-radius rela- 
tions, depending on three functions describing the effect of rotation on star structure. These 
relations are exact in the sense that they have been obtained without any special assump- 
tions. By assuming some appropriate forms for the unknown functions we have obtained a 
general description of the mass and radius of the rotating neutron or strange stars, which 
generally and for a broad class of equations of state can reproduce the values obtained by nu- 
merical integration of the gravitational field equations with a mean error of around 5%. The 
expressions of the unknown parameters have been chosen following a close analogy with 
the static case,whose relevance for the study of rotating general relativistic configurations 
seems to be more important than previously believed. These functions also incorporate some 
other general relativistic effects not explicitly taken into account, like the variation of the 
moment of inertia of the star with the angular velocity. Ravenhall and Pethick [19] have 
presented a formula valid for a broad range of realistic equations of state of dense matter 
expressing the moment of inertia in terms of stellar mass and radius. We have not used 



28 



these results,obtained in the slow rotation limit because, at least in the case of strange stars, 
the Newtonian expression of the moment of inertia also leads to a quite accurate physical 
description of the rotating objects. As an application of the obtained formulae we had given 
a derivation of the "empirical" formula relating the maximum angular velocity to the mass 
and radius of the static maximal stellar configurations. 

The possibility of obtaining the basic parameters of general relativistic rotating objects 
in terms of static parameters could lead to a major computational simplification in the study 
of rotation. The relation between the presented formalism and the Einstein gravitational field 
equations will be the subject of a future publication. 
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FIGURES 




FIG. 1. Variation of the radius (in units of 10 6 cm) of the static strange star as a function 
of the parameter r\. The solid curve represents the values of the radius calculated with equation 
(|7(i|) , while the points represent the values of the radius of the quark star obtained by numerically 
integrating the general relativistic structure equations of the star. 
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FIG. 2. Variation of the mass of the static strange star (in solar mass units) as a function of the 
parameter 77. The solid curve represents the values of the mass calculated with equation ([FT]), while 
the points represent the values of the mass of the strange star obtained by numerically integrating 
the general relativistic structure equations of the star. 
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FIG. 3. Variation of the total energy (matter+gravitational) of the static strange quark star. 
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FIG. 4. The radius of rotating strange star (in units of 10 6 cm) ,as a function of the angular 
velocity, for values of the central density in the range r] G [5, 16.63] .Values obtained from the 
analytical expression (solid curve) are compared with data obtained by Stergioulas et al. [42] 
numerically integrating the gravitational field equations (points). 
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FIG. 5. Variation of the radius of the rotating strange star (in units of 10 6 cm) as a function of 
the angular velocity f2(s -1 and of the central density rj. 
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FIG. 6. The mass of rotating strange star,expressed in solar mass units, as a function of the 
angular velocity, for values of the central density in the range r\ € [5, 16.63] .Values obtained from 
the analytical expression (solid curve) are compared with data obtained by Stergioulas et al. [42] 
numerically integrating the gravitational field equations (points). 
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FIG. 7. Variation of the mass of the rotating strange star, expressed in solar mass units, as a 
function of the angular velocity and of the central density r\. 
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